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Abstract 

We develop a method to obtain the general solution of the Laplace 
equation in d-dimension in ultraspherical coordinates. 

The coordinates of a point in d-dimensional space are given by the position 
vector Yd = {xd ), |r<i| = r. We define the ultraspherical coordinates 

by Tj = rj{cos6j\i.j + sniOjVj^i) = x^kj + rj-i, < 6'j < tt, j = 3, • • • , d, 
1*2 = T2 cos0ki + r2 sin0k2, < (/> < 27r, k^ • kj — 5ij. The hat symbol over a 
bold letter always means a unitary vector. 

We want solutions that are finite in the angles and the frontier conditions 
are radial only and the domain in the angles are complete, i.e, < (/) < 27r and 

< 0j < TT, J = 3, • • • , d. 

1 The laplacian operator 

Before to find the solution of the Laplace equation in ultraspherical coordinates, 
we want to find the explicit form of the laplacian operator. Consider the operator 

dz^ ^ p" dp \ dp) ~ dz^^ dp^^ pdp' ^ ' ^ 

Applying a coordinate transform z — r cos 6 and p — r sin 9: 

d „ d cos 9 d 
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For d = 2: 

^' = d^ + d^ = -pd-p[%)+7W ^ = ^'^°^'^' y = (1-3) 

For d^3: 

- — V^-— -—( —\ — — 
Using (II. 2p with a = 1 and z = rcos9 and p = rsin0: 

(1-4) 

By induction, supposing valid for d — 1: 

vLi = ^^'-i7f^) + ^^'^-i- (1-5) 

Using (II. 2p with a — d~ 2 and = r cos^^i ''rf-i = 7'sin0(j we have 



were 



' (1.6) 



1 5 / . . 2. a \ . 1 



sm" 25/^ 96*^ V dOdJ sm^Od 

We see from p.6p that we can obtain the solution of the Laplace equation if we 
known the eigenfunctions of the angular operator Ld- 

2 The general solution 

Let the eigenfunctions of Ld given by ^'(ri), where is the angles coordinates 
n={9d,0d^i,---,O,<P). 

Ld-^pin)^ (3-9pin). (2.1) 

Now supposing the solution of the Laplace equation by the separable form 
^p{r, = R/3{r)^fj{Q), we get the radial equation: 



r'^ ^ dr \ dr J 

This is a second order Euler equation whose solution is given by r", where a is 
to be determined. We find 

(3 ^ ^a{a + d~2). (2.3) 
Then the general solution is the form 

^(r, n) = ^(^(a)r" + B{a)r-^"+'^-^^)-^a.{n), (2.4) 

a 

where L<j*Q(il) = -a{a + d - 2)\E'a(ri). 
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3 The eigenfunctions of the angular operator 



Since we are interested in the general solutions that the frontier conditions 
are radial only, the angular part are the same for all type of frontier condition. 
Consider first the non-trival solution of (j2.4p which is independent of the angles. 
Clearly a solution is for a = with A{0) — 0, that is B{0)r^^'^~^\ To construct 
a solution with an angular dependency, we observe that the laplacian operator 
is invariant under translation. Making a unitary translation in the Xd direction, 
we conclude that 



\rd-kd\'^-^ (l + r2-2rcos0d)'^ 

is a solution of the Laplace equation for r ^ 1 that goes to zero when r — > oo. 
We define the ultraspherical polynomials by 

oo 

Y,r'PiAcos9d), r<l. (3.1) 



(1 + — 2r cos 6a 



d-2 



Note that for d = 3 the equation above gives the definition of the Legendre 
polynomials. 

Applying the laplacian operator (|1.6p in both sides of (|3.ip and taking the 
fact that r' are linearly independent for each I we have 

LdPi4{cos9d) = + d- 2)P,,d(cos0rf). (3.2) 

For d = 3 the equation above reduces to the Legendre equation. Note that this 
eigenfunctions depend only of the angle Od- This solutions can only be applied 
in problems that are independent of the other angles. 

To find a general solution for d-dimension, observe that a solution in {d— 1)- 
dimension is also a solution in c?-dimension. Consider first the case d = 3 for 
a better understand. Is well known that — p^e*™"/" is a solution for d — 2. 
Then is also a solution for d = 3. But linear combinations of this solutions can 
solve a part of two dimensional problems only (finites for p = 0). To find a 
tri-dimensional solution let us try 

where a is to be determined by the requirement that fm being a solution of 
the Laplace equation. Note that the dependency of tri-dimensional angle in the 
denominator and for m = and a = I we have a solution of the type (|3.ip with 
d = 3. Applying the laplacian operator in both side of (|3.3p we have 



since V3 — + V2 and r — = r2 + (.T3 — l)k3, with p = r2 = rsint 
X3 — r cos 6*, we get 



72 1 ^ «("~ 1) 

^|r-k,,h |r-k.,|"+2- 



Thus we have 



'31 

^lf^ = -a{2m+l-a)y-^^^ = 0. (3.4) 
There is a non-trivial solution for a 7^ in (|3.4p . that is a = 2m + 1. Then 
_ r^e"""^ _ r" sin™6'e™'^ 

J"^ ^ Tz r. i2m+i ^ 77" ^ ;:; TTiii+r' 

|r — K3I ^ (1 + — 2r cost') 2 

is a solution of the Laplace equation in d = 3 dimension. To obtain separable 
solutions, let us take the m derivatives with respect to x = cos0d in p.ip : 

, f -±r'i^. ,3.) 

(1 + r2 - 2rx) —2 — ^ da;™ 

where a(m, d) = (d — 2)(i((i + 2) • • • {d + 2m — A). Taking d = 3 in p.6p we obtain 

r-E'^ ' x = cos0,r<l. (3.7) 



2m + l 



(l + r2-2ra;) - ;^(, 
Substituting p.7p in (|3.5p we conclude that 

V r' sin"0— -ipe™^ rr = cos^, r < 1, Pi = P,,3 (3.8) 
1=0 

is a solution of the Laplace equation in d — 3. Let us define the associated 
Legendre function by 

P;'"(cos6i) :^sin"6> ' , x^cosB. (3.9) 

Note that since P;(cos6') are polynomials of degree Z, we have that < m < Z. 

By applying the laplacian operator in the equation (j3.8p and taking the 
account the linear independence of we have 

h (Pr(cos0)e™^) = -l{l + l)Pricose)e'""f'. (3.10) 
Writing L3 in the explicit form, wc have the equation satisfied by the P™ (cos 6) : 

1 "^.ta.^rVff^ + D-^ilpr^o^ <Mr> 



sinOde \ dO 



This is the associated Legendre equation. It is a Sturn-Liouville equation and 
this functions P™ are orthogonal in < 6* < tt with a weight sinfl for a fixed m: 

( sin 9 Pi"" PiT d9 ^0, I ^ I' . (3.12) 
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Since that e""* are orthogonal in [0, 27r] we have a basis of functions that are 
orthogonal in [0, 27r] x [0, tt]. We have 21 + 1 linearly indepependcnt functions 
for a given /, since P™(cos 0)e~*'"'^ also satisfies p.lOp . 

In order to obtain the general eigenfunctions of Ld, let us suppose by hy- 
pothesis that a solution of the Laplace equation in {d — l)-dimension is given 

by 



W2-™i =Cr*™.-2-™i(^d-i---,0,</>), 03 = O,92=(I}, (3.13) 



where 5'm<i_2 - mi (^d-i • • • , ^, is a function of the angles of (d — l)-dimension 
only. For example for d = 3 a solution in 3 — 1 = 2 is the form g„i — r^e"""^. 
Let us try a solution like (|3.3p : 

Jma-2---m,i — s- . ■ {•^■^^J 

|r-kd|« 

By the hypothesis 

^ d9md-2---mi = ^ d—l9rnd-2---mi — 0, 
Td-l ■ ^ d-l9md-2---ini — ™d-2.gmd-2 ■ ■ -mi , 

„2 . - n ^ a{a - 2md-2 - d + 2) _ 

|r-k<i|"+2 

For a 7^ we have a non-trivial solution with a = 2md-2 + d — 2: 

/md_2--mi - 1^ _ ^^|2„i^_2+d-2 ~ | j. _ j^^ 1 2md-2+rf-2 ^' 

Now using p.6p with to = md-2 we conclude that 

^ sin"-^0,^-;;^^^M'™,_,...^, • • • , 0, 0), X - COS r < 1, 

i=o 

(3.16) 

is a solution of the Laplace equation in d-dimension. 

We define like (|3.9p the ultraspherical associated Legendre function by 

d™- Pj j(t^ 

Pl^^{cosdd)^sm"'dd^^, x = cos0d. (3.17) 

Applying recursively for d — 2 until d = 3, we have the form of the angular 
function ^P: 

^lmd.2-m2rr., = P™'"' (^OS ed)Fr;-^_ i (cOS i ) • • • P^^COS 0)e™^^ 

(3.18) 

where < TOd-2 < ^ < to^-s < rnd~2, ■ ■ ■ < mi < TO2. 

Taking the laplacian operator in p.l6p and the linearly independence of r' 
we have 

-^d^lmd-2---'m2mi — l{l + d 2)\I/^yj^^_2...rM2rMi ■ (3.19) 
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To obtain the equation satisfied by the P[^{cos9d), we make nid-2 = ™ and 
mi = 7712 = ■ ■ ■ 77id-3 — and use the exphcit form of Ld'- 

1 d O-2^,^Uf;(^ + rf_2)-!!^(!!i±^V."=0. 



^^T^ddO-dV '"'^J^K' ' sin^^, 

(3.20) 

Note that for d = 3 this equation reduces to equation of the associated Legendre 
functions. This equation is a Sturn-Liouville type and the P;™ are orthogonal 
with weig ht sin'^-^Od in the interval [0, tt] for a fixed m: 

f sin''-^edPrdicosed)PiTdicosed)ded^o, i^i'. (3.21) 

Jo 

Now using p.lSp we have the orthogonality for the eigenfunctions '^imd-2---rn2mi ■ 



(3.22) 

where dil = sin'^^^^d sin'^^^^d-i • • • svaOddddOd-i ■ ■ ■ d9d<j). 

We have a set of mutually orthogonal functions of entire domain of the 
angles. The number of this functions for a given I is 

_ {d + 2l^2){d + l-3)l 
^'^^^ {d^\ ■ ^^-^^^ 

4 The Laplace equation solution with radial fron- 
tier condition 

Using the results of the previous sections we conclude that the general solution 
of the Laplace equation in d-dimension with radial frontier condition and finite 
for all value of the angles is 

$(r,r!) = ^ E •••E E ^'™---'n.W^;™._....™,(l^)+C.C, (4.1) 

l—0mfi_2—0 m2— Omi— 

where 

F, (r) - Ai r'- + R, -{i+d-2) 9) 

^ lmd~2---^l\ / — i^d — 2 ■■■mi ' ^ -'-*lmd-2^ ■ ■mi ' i v / 

- ^imd_2 - mi are complex constants to be determined by the frontier 
conditions and C.C is the complex conjugated. 

To write (|4.ip in a compact form, let us define the ultraspherical harmonics: 

^^id-2"""^l — flmd-2^ ■ ■mi^ lmd-2^ ■ ■mi : (^■*^) 
Ylmd-2 mi — /imd_2---mi^;md_2---mi' (4-4) 

where /im<i_2 - mi is a normalization coefficient such that 



n 



■^')n^_2---m'i ~ '^''"^"id-2"i^_2 ' ' ' ^m2m'^Smim[ i (4-5) 

and are given by 

f, - /v('^) ivid-^) .../V(3) (AQ) 

./imd_2---mi — ^*;md-2 nid-2md-3 ^^m2mi^ \^-") 
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where -/V,lm is the normaUzation factor of the ultraspherical associated Legendre 
function in d-dimension: 

{Ni^J^y r sin'^-'e, (P„™ (COS0,))' d9, = 1, (4.7) 



(,)^ l {2n + d-2) a^^i (rf-3)!(n-m)! 

"™ Y {d-2) ild {d + n + m-3y: ^ ^ 



Now (|4.ip can be written as 



ms m2 



<P{r,n)^J2 E ■■•E E (4-9) 



=0 md-2=0 m2=0mi= — m2 
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A The volume element in ultraspherical coordi- 
nates 

We see from the definition of the ultraspherical coordinates that 

Xj — Tj cos9j, Tj-i = rj sin6'j, < Oj < tt, j > 3, (A.l) 
x = pcos(j), y = psincf), < (/) < 27r, p — r2- (A. 2) 

The volume element in d = 2 is dV2 = dxdy = pdpd(f). In d = 3 we can write 

dVa = dx^dV2 — pdx^dpd^i. (A. 3) 

Now using (|A.ip with j = 3, .13 = r^cosd^, p — r2 — r^sind^ we have that 
dx^dp — r^dr^dO^. Then 

dVz — (ra %YiiB-^ri,drj,dQz<^4> — t\ sin O^dr^dOsdc/) ~ r^dr^dils 

Now supposing that the volume element in {d — l)-dimension is 

dVd-i ^ r^aZldrd-idfld-u (A.4) 

we have 

dVd = dxddVd-i = r'^Zldxddrd-idVld-i. (A. 5) 

Now using (IA.1|) with j = d, Xd — Vd cos Od, rd^i — Vd cos 9d, we get dxddvd-i = 
TddrddOd- Hence 

dVd = r^^i sin'^-^eddrddeddnd-i = ^-^drddnd, (A.6) 
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where 

dUd = sin'^-^0dd9ddnd^i. (A.7) 
Using (jA.7[) recursively, we conclude that 

dfld = sm'^-^Od sin'^^^e'd-i • • • smOsdedded-i ■ ■ ■ dOadcj). (A.8) 

The solid angle ild is given by 

nd=( f sin^-^eddOd) nd-i = V^^^^d-i = 2-^. (a.9) 
B The normalization factor of P'l^ 

From the definition of the ultraspherical polynomials we have 

1 °° 

Ir _ ^ E ;i+^^M(cos7d), (B.l) 

where cos7ti = f • f ' — cos 6 d cos 9'^ + sin 9d sin 6'^ cos jd-i and r'<(r>) is the 
smaller (larger) of r and r' . Taking the laplacian in both sides of (|B.1|) we 
obtain 

(2 - dW^^ir - r') = £ (0.(r)^ - ^^^^t^^) PiAcosjd), 

(B.2) 

where S'^{r — r') is the Dirac delta distribution in d-dimension, 

r'^ ^ sm" ''Bd smOs 
and Od{r) is the radial operator, 

> - rd-i dr y dr 

Multiplying (|B.2p by r'^~^ and integrating from r' — e to — e with e ^ 0, we 
have 

(d- 2)f^,fc^ . . . = f;(2/ + d- 2)P,,(cos7.). (B.3) 

sm" ^Od sm6'3 ^ 

Since the ultraspherical harmonics are orthogonal in the domain of all angles 
we can write 

I m2 

Pi.d (cos 7d)= ^ ••• ^ Ai Yl„^,_2-mA^d)- (B.4) 

md_2— mi— — m2 

Using the orthogonality condition (|4.5p and (jB.3p with (jB.4p wc obtain 

_ (rf-2)», 

^<md_2---mi — 2^ _|_ _ 2 "''^^""'"^ 
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Thus we obtain the ultraspherical harmonics addition: 

(B.5) 

We can write the equation above in a useful formula 

PiA^os-ta) = Ki^d E (2mrf_2+rf-3)(Af/„t_,)'Pr/"'(cos0d)P,7-^(cos0:,)P,„,_,^,_i(cos7d-i), 

nid-2=0 

(B.6) 

where 

^ rid (t^ - 2) 

Now deriving the equation (IB.6|) with respect to 9'^ n-times with 6*^ = and 
equalling the independent terms we have 

^(2n + d-3)X,,P,(;)(l)P„5^,(l) 

where Pj|;"^(l) and Pn\i-i{^) is the n-th derivative of Pi^d{x) and P„.(i_i(x) 
with respect to x evaluated in x = 1. This can be obtained from (j3.6p with 

" ^' '(„) _ [d - 2)d{d + 2) ■ ■ ■ (d + 2n ~ 4)(rf + n + / - 3)! 

'■■^ " G-n)!(d + 2n-3)! ' ^ ' 

Then 



p'{")n^p'(«) (d + n + /-3)! 

^'•^ ^ ' "■''-1^ ^ (d-4)!(?-n)!(d + 2n-3)' 



We finally obtain 



(rf)^ {21 + d- 2) (d-3)!a-n)! 

V (d-2) rid (d + ? + 7i-3)! ^ ■ ' 
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